A SUPEREXPONENTIALLY CONVERGENT FUNCTIONAL-DISCRETE 
METHOD FOR SOLVING THE CAUCHY PROBLEM FOR SYSTEMS OF 
ORDINARY DIFFERENTIAL EQUATIONS 



Makarov V.LQ and Dragunov D.V]^] 

Department of Numerical Mathematics, 

I Institute of Mathematics of NAS of Ukraine, 

3 Tereshchenkivs'ka Str., Kyiv-4, 01601, Ukraine 

o 

<L> MSC: 65L05, 65L20, 65L80. 

Q 

In the paper a new numerical-analytical method for solving the Cauchy problem for sys- 
tems of ordinary differential equations of special form is presented. The method is based 
i— i on the idea of the FD-method for solving the operator equations of general form, which was 

proposed by V.L. Makarov. The sufficient conditions for the method converges with a su- 
perexponential convergence rate were obtained. We have generalized the known statement 
pH about the local properties of Adomian polynomials for scalar functions on the operator case. 

Using the numerical examples we make the comparison between the proposed method and 
the Adomian Decomposition Method. 

Introduction. A great number of papers published during last two decades are de- 
voted to the Adomian Decomposition Method (ADM) and its applications (see pQ [5] 
and the references therein). This method was proposed in the middle of 80-th by Amer- 
ican physicist George Adomian. It has aroused great interest among mathematicians all 
over the world. The very foundations of this interest lie in the fact that this method 
belongs to the class of analytical methods: it allows one to compute approximate solution 
of nonlinear operator equation by solving a sequence of linear equations. Using the ADM 
we looking for the approximate solution in the form of partial sum of the series which 
converges to the exact solution. Due to its analytical property the ADM allows one to ob- 
tain the approximate solution even in the case when initial operator equation depends on 
the parameters. In other words, the ADM allows us to approximate the multi-parametric 
family of solutions in the functional form. 

The general idea of the functional-discrete method (or simply FD-method) is quite 
similar to the ADM's one. The FD-method was for the first time proposed in [6] where 
it was applied to the Sturm-Liouville eigenvalue problem. In [6] it was showed that 
the convergence rate of the method increases as far as increases the order number of the 
eigenvalue which we are going to approximate. So, the FD-method demonstrated amazing 
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results that are unachievable for the finite elements method, mesh method and other 
purely discrete methods. Leter, in [7j [11] , the FD-method was applied to the nonlinear 
eigenvalue problems and the sufficient conditions that guarantee its superexponential 
convergence rate were found. Applying the FD-method to the boundary value problem, 
the authors of [12] proved, that in this case the superexponential convergence rate of the 
method can be reached (under certain sufficient conditions) as well. In addition to this, 
in [T2] the general scheme of the FD-method for some class of operator equations in an 
abstract normed spaces was presented. In [21] the FD-method was applied to the Cauchy 
problem for ordinary differential equation of the first order and the sufficient conditions 
which guarantee superexponential convergence rate of the method were found. 

Both the ADM and the FD-method are originated from the homotopy method. Both 
this methods approximate the exact solution of the nonlinear problem using the partial 
sums, which terms are to be found as the solutions of some sequence of linear problems^} 
But the essential difference between those methods is expressed by the fact that the FD- 
method has a built-in parameter — the parameter of discretization. By varying this 
parameter we cane achieve the convergence of the FD-method even in the cases when the 
ADM is turn out to be divergent. The last fact will be demonstrated below by several 
numerical examples. 

In the present paper we offer the algorithm of the FD-method for solving the Cauchy 
problem for systems of ordinary differential equations (SODE), which are based on the 
general scheme of the FD-method for operator equations [T2]. It is the further development 
and generalization of the theoretical results obtained in [15] . The main result of the paper 
is presented in theorem |2j which contains the sufficient conditions for the FD-method 
applied to SODE converges with superexponential rate. 

Let us consider the following Cauchy problem for SODE of the first order 

^lt(t)-N(t,lt(t))Tt(t) = ~$(t), lt(t ) = lt eV m {R), te[t„,+oo), (1) 

where N (t, it) e M m (c°f ([t 0) +oo) x R™)) , it (t) = [ Ul (t) , . . . , u m (t)f , ~$ (t) G 
V m (C ([t , +oo))) . Here and below we use the notation M m (F) to describe a linear space 
of square matrices of the m-th order with elements belonging to a linear space F. Similarly 
to this, by V m (F) we denote a linear space of column vectors of the m-th order with 
elements from F. We also assume that problem ([I]) possess a unique solution. 

Using the FD-method for solving the Cauchy problem ([!]), we looking for the approx- 
imation of the exact solution u it) in the form of the partial sum (see [T2"] ) 

lt(t) = J2^ {t) (t), P eN\J{0} (2) 

3 Generally speaking, this is not true for the FD-method, because the base problem is not linear one (see later). 



2 



of the convergent series 



(3) 



i=0 



In this case the nonnegative integer number p is called the rank of the FD-method. To 
emphasize this detail we usually say: "the FD-method of the p-th rank" . 

To apply the FD-method to the Cauchy problem ([T| we need to introduce a grid on 
the interval [to, +00) : 



Q = { t < h < h < ■ 



lim 



t n = +00 \, h = sup {hi = ti- ti-i} . (4) 



After that we can consider the following recursive system of linear Cauchy problems with 



respect to the unknown terms if) = up (t) Um if) , j G N |J {0} of the series 



d (t) - N (t, lt<® (*<_!)) ^(°) (t) =~${t), te U) , Vi G N, 



dt 



(5) 



^(o) (to) = to, P (0) (t)] t=u = {U + 0) - ^(°) ft — 0) = if, Vi G N. 



^) (t) - N ft l^W ft_0) ^W) (t) 



dx 



771 X 



^(°) (t) + (t) , te [t,-i, u) 



(6) 



^0) ft) = [it® (t 



t=u 



Vz,jGN, 



where 



FW ft = £ V P (N ft (■)) ; ^ (0) ft-i) , • • • , ft-i)) ft + 



p=l 
J'-l 



+ E I^'-i-p (n ft (■)) ; ^ (0) ft , • • • , ^°'- 1 - p) ft) - 

p=0 

-Aj-i- P (N ft (■)) ; ^ (0) , • • • , ft_x))] ft + 

(N ft (•)) ; ^ ft_i) , . . ■ , ft_ x ) , <°) (t) , 



(7) 



te[ti-i,ti), Vz,jgN. 



Here A k (N ft (•)) ; ^0, ^1, • • ■ , ^fc, ) = ( t, £ r^, 

i=0 



denotes the Adomian's 



T = 



polynomial of the k G N|J {0} order for the operator N ft (•)) ( see [1], [3]). 
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Problem ^ is called the base problem. It can be considered as an example of the 
Cauchy problem with piecewise constant argument (see [13J). If we denote by Ui (t) = 
Ui (t, N (t, ^(°) (ti-i))) e M m (C 1 [U-uU]) , i G N the solution of the following Cauchy 
problem 



|^ (t) - N (t, l^W ^ (t) = 0, Ut (t<_i) = t6 [<<_i, t*] , (8) 

the solution of problem ([5]) can be expressed in the form 

t 

(t) = U (t) (vO + y A, (t, (0 de, t e t«] , Vi e N, (9) 

t<-i 

where 

A, (t, = ^ (*) u; 1 (0 , te ti] , Vi e N, (10) 

is the Cauchy matrix (see, for example, [16, p. 412]). Similarly we can express the 
solutions of problems (pi), (m) Vj G N: 



C0 (t) = 17, (t) W + f K (t, T< (0 (**_!) d£+ 



+ f K(t, F<a (0 de, * e *i] , Vi g n, 



where Tj (£) G M m (C 1 denotes the following matrix 

(t) = 



(? N (t, ^ (ti_x)) tf<°> (t) , . . . , (t, tf<°) (t,_0) ^ (0) (t) 



■ i v'" " vf— i./,' " v / -i - : - " v " v' — •• v / • \ / 

oui aw,; 
Vi G From expressions ([9j, (11) and (]7| we can easily obtain the following 

statement. 

Proposition 1. TTie column vectors (t) , (t) , . . . , which are the solutions of the 
recursive linear Cauchy problems (|5])-(|7j), always exist and are unique on [t , +oo] . 

Definition 1. We will say that the FD-method for the Cauchy problem (fij) converges (to 
the exact solution of problem Q) on [to, to + H) , to < H < +oo, if there exists a positive 
constant h G K suc/i £/iat /or every grid a) Q satisfying h < h, series with terms 
obtained as the solutions of problems (Jsj) — (J^J) ^ converges absolutely and uniformly (to the 
exact solution of problem ([B) on [to, to + H) . 

The local properties of the Adomian's polynomials. Let Ei and E 2 be some 
normed linear spaces over R with norms ||-||^ and ||-||^ respectively. A set of linear 
operators which map Ei to (on) E 2 we will denote by L (Ei, E 2 ) . Lets consider an operator 
N (u) (nonlinear, in general), which maps a given open subset G of the space Ei into E 2 , 
i.e. N: Q — > E 2 . Let us recall some definitions of important concepts, which we will need 
later ( see [M p. 28-30]). 
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Definition 2. An operator N"( fc ) (ui, U2, . . . , u k ) : E^ — > E2 is said to be a k-linear op- 
erator, if it is linear with respect to each argument Uj G Ei, % G 1, k. 

An operator (u 1; u 2 , . . . , u k ) is called to be symmetric, if its value is independent 
on the order of its arguments. 

We will denote the set of /c-linear operators (ui, U2, . . . , u^) : E^ — > E2 by 

A (Ei,e 2 ) . 

Definition 3. A given k-linear operator (ui, u 2 , u k ) is said to be bounded, if 

there exists a constant M such that for every Uj G Ei the following inequality holds 

||N( fc ) (u 1; u 2 , . . . ,u fc )|| (2) < M || Ul || (1) ||u 2 || (1) . . . ||u fc || (1) . (13) 

For the smallest constant M which satisfies inequality (13) we will use notation ||N( fc )|| . 
It is called the norm of the k-linear operator N^ fc \ In other words, 

||N«||= sup ||N( fe )(u 1 ,u 2 ,... ) u fc )|| (2) . (14) 

||ui|| (1) <l,-,||u fe ||«<l 

Definition 4. An operator N (u) is said to be differentiable by the Frechet at the 

point u G G if there exists a bounded linear operator (u ) G C (E l5 E 2 ) such that 

N (u + h) - N (u ) = N« (u ) h + u (u , h) , (15) 

where ||w(uo,h)|| (2) = o(||h|| {1) ) , ||h|| (1) 0, Vh: Uo + h G G. The linear operator 
N^ 1 ) (u ) , when it exists, is called the 1-st Frechet derivative of the operator N (u) at 
the point uq. 

Using the induction we can define the k-th Frechet derivative of the operator N (u) at 
the point u , Vfc 6 N [T71 p. 262]. Namely, if the 1-st Frechet derivative of the operator 
N (u) exists at each point of the certain open neighbourhood Gi C G of the point u , 
then we can consider the mapping (or the operator) 

N (1) (u) : Gi ->£(Ei,E 2 ). (16) 



If there exists the 1-st Frechet derivative of operator (16) at the point u , then we will 
denote it by N^ 2 ^ (uo) and call it the second Frechet derivative of the operator N (u) 
at the point u . It is obvious that (u ) G £2(E 2 ,E 2 ). In general case, the k-th 
Frechet derivative of the operator N (u) at the point u belongs to the space (Ej, E 2 ) . 
It is denoted by N^ fc ^ (uo) and is defined as the 1-st Frechet derivative of the operator 
N^" 1 ) (u) : G fe _! -± (E^ X ,E 2 ) at the point u G G k ^. 

Definition 5. An operator^ (u) is said to be k (k G N) times differentiable by the Frechet 
(or in the Frechet 's sense) on the certain open set G k C G, if there exists the k-th Frechet 
derivative of the operator N (u) at each point of G k - 
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Below we give the definition of the Adomian polynomial of the n-th order for the 



A n (N(.);{u l ]l Q )=-—N(<t>(T)) 



(17) 

r=0 



operator N (u) with respect to the "variables" u, G E, i G l,n. This definition is more 
detailed then the definition from jH p. 59]. 

Definition 6. Let the operator N (u) : G — > E 2 is differentiable by the Frechet at the point 
uq G G up to the n-th order inclusively, n G N. Additionally, let $ (r) = f $ (r, Ui, . . . , u n ) 
be an operator which maps K x E" (r G R, Uj G Ex, i G 1, n) to Ei and is defined by virtue 

n 

of the formula $ (t) = E t 1 U{. Then the Adomian polynomial of the n-th order for the 

operator N (u) is defined to be an operator A n (N (•) ; [uj]™ =0 ) = A n (N (•) ; u , Uj, . . . , u n ) , 
which maps E™ to E 2 and can be explicitly obtained by the formula 

1 d 
nldr 

If n = then we assume 

A (N(-);u ) = N(u ). 

For convenience we will use the following notation for the fc-linear operator (uo) G 
Ck (E*,E 2 ) 

(fcl fc2 
ui, . . . ,ui,u 2 , .7. ,u 2 ' " v "" 1 

^#)(u )(kk... 1 n) 1 (18) 

i 

Vui G Ei, h G N (J {0} , i G 1, 1, k j = k - 

Theorem 1. (see [U p. 60]) Let an operator N (u) is differentiable by the Frechet at the 
point uq G G up to the n-th order inclusively, n G N, then Vuj G Ei, Vi G l,n 




1 /pi 
A„(N(.);[u^ =0 )= £ - T — T N^ + -+^ (uo) [ui] 

/>1 Ai- \ 



Pn 



(19) 



2 kp k =n 

p fc eN|j{o} 



1 [ Ot\—Ot-2, «n-«n+l 



(20) 



A„ (N (■) ; [u,^) = £ ii N(«0 (u ) ( [uj , 

;^on(«.-%i)! 

Lets consider the subspace £i of the normed space E 1; which is also a normed space 
with the norm ||-|| . Let ||-||i be a scalar function defined on Si, which satisfying the 
conditions of seminorm (see, for example, p2J p. 564]), then a scalar function |||-||| , defined 
by the equality 

|||u||| = max|||u|| (1) , llull^} , Vu G S u (21) 
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is a norm (see, for example, [T71 p. 41-42]) in the space S\. 

The following lemma is a generalization of lemma 2.1 from [T2] . 

Lemma 1. Let the following conditions are satisfied 

1 ) An operator N (u) is n times (n G N) differentiable by the Frechet on the open 
convex subspace G n of the space Ei, G n f]Si ^ and there exists a scalar function 
N (u) G C n (R) which satisfies the following inequalities 



|N«<u)|| 



VfcG0,n, Vu G G n . 



(22) 



l(!) 



2) For some positive constant h G R and operator P, : £ i — > Ei i/ie following conditions 
hold true 

||u-Pfc(u)|| (1) </i||u||«, Vuefi, 
P/,(u)eG B , VueG n f|£i, 
||P, (u) - 0! (P h (u) - u)|| (1) > ||P, (u) - 9 2 (P h (u) - u)|| (1) , 



(23) 
(24) 



(25) 



V0i,0 2 eM: 0i > 2 > 0, Vug£i. 



Then 



A k N (•);[< 



(2) 



<A fc iV(-); 



U; 



j=0 



(26) 



Vu G G n , Vuj G Ex, iel,k, Vfce0,n-1, 
Ak (N (■) ; [u,]to) - A k (N (•) ; [P h (uA}^ ® 

<hA k (n® (•) x (•) ; [IKlto) , 



< 



(27) 



Vu GG n U£i, VuiG^i, ie l,fc, V£;G0,n-2. 
Proof. Let the conditions of the lemma are fulfilled. At a first step we are going to 
prove inequalities (26). Lets fix an arbitrary k G 0,n — 1. Using theorem [TJ specifically, 



representation (20), as well as inequalities (22), we would obtain 

llN( tt fWlll 

(2) 



A k (N (•) ; [uiJto 



||N(-)(u )||n(||u 



< 



E 

ai+...+a fc =fc 
Qi>...>a fc+1 =0 
« 8 6N|J{0} 



i=l 



< 



n («i - a i+ i)\ 

i=l 



i\r( ai ) 



< 



E 

ai+...+a fc =fc 
ai>...>a fc+ i=0 
ai6N|J{0} 



U 



n (k 

i=l 



Cti— C«i+1 



(28) 



n («i - ai+i' 
i=i 
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U; 



i=0 



Vu G G n , Vuj G E 1; i 6 1,1c. 



Let us prove inequalities (27). For this reason we fix an arbitrary k G 0, n — 2 again. 



Using representation ( 20 ) and the generalization of the Mean value theorem on the oper- 



ator case (see, for example, [T7, p. 248]), as well as the properties of the operator's N (u) 
derivatives, we will obtain: 



A k (N (•) ; [u.to) - A k (N (•) ; [P, ( Ui )t 



E 

a 1 +...+a k =k 
ai>...>a fc+1 =0 
"i6N|J{0} 



'a.\ — a.2 a k~ oth+i 

N^>(u ) ( M [u fc ] 



I] (oci - a i+1 )\ 

i=l 

. _ '-k+l " 

N^)( Ph ( Uo) ) [p ft ( Ul )],...,[P h (u*)] 



«1- C*2 Q fe — Q fe+1 



(29) 



E 



ai>...>a fc+1 =0 II («i - «i+lj! 
ai6N|J{0} i=l 



Q1-CI2 Q fc— "fc + 1 



1] (atj - a i+ i)\ 
i=i 

{Nt«) (uo) - N<«0 (P h (u ))} (JuJ, . . . , "[u k )"j + 

, / ax — a% — % i— 1 «2— "3 a fc — a fe+i\ 

^ N^P^uq)) (ux-P^m), [ Ul ] ,[P,( Ul )], [u 2 ] [u fc ] J + 

+ ]T N^(P,(u )) [P,( Ul )],u 2 -P,(u 2 ), [u 2 ] ,[P,(u 2 )],..., [u fc ] + 



i=l 



ax— a.2 



. . , a^—ii—l 

. . . + J2^ (ai) (Ph K)) ( [Ph (Ui)], . . . , [P fc (Ufc-0], U fe - P h ( Ufc ) , [U fc ], [P h (Ufc)] 
i=l 



From condition (22), the convexity of the set G n and from properties (23) - (25) of the 
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operator P h (u) we would obtain 



ai-a 2 "fc-Ofc+l 

{N(«0(uo)-N<«)(P fc (uo))}( [m],..., [u k ] 



|N^ +1 ) (P h (u )-e(P h (u )-Uo))| 



(2) 



xiKuo-p^uo^ii^nfiKuwr^ 1 -. 



i=l 



< /iiV^+ 1 ) ( ||P h (u ) - 9 (P h (u ) - u ) 



luolfniK 

i=i 



iW 



< 



< Ml Uo 



K 1 ) 



luoii^nniui 



< 



< hN^ (llluoll) luoll EI (|ui|) Q! - Ql+1 , G (0, 1) . 

i=i 



Taking into account (29) and (30), we would get the inequality 



A k (N(-);[u,t )-^ (N (•) ; [P h (u 4 )t ) 



(i) 



< 



1 

^ E Hi 

Ql + ...+«fc=fc TT / \ I 

"l>--->«fc+i=0 11 [ a i ~ Oii+l)- 
Q<£NU{0} i=l 



WV^+D (luoll) llluoll n (INI 
1=1 



+ 



|N^> (P^uo^Hllux-PfcWII^n (IK 



i=2 



K 1 ) 



X 



E fl|ui|| (1) 



1=1 



(llPft(ui) 



1(1) 



t-1 



fc-1 



+ ||N^) (P, (u ))|| \\u k - P h (u fc )|| (1) n ( IIP/. (Ui 

i=l 



xlEf||u fe || (1 T fc "(H p /'("'i 



U=l 



l(l) 



i-1 



1(1) 



< 



C*i— Oi + l 



X 



i 



n (K 



X 



Ofl + ... + Qfc=fc -| — r / \. 

^i>...>a fc+ i=0 11 [&i — a> i+ i)l 
a,eN|J{0} i=X 



i=l 



X 



jy(ai+l) 



+ «iiV( ai ) 



hA k ( JVW (•) x (•) ; [|||u 



* III J i=0 



Vu G G„ U Ex, Vui G fi, i G 1, fc, 
which we needed to prove. The theorem is proved. ■ 



To prove the main result of the paper we need to use the partial case of lemma [T] which 
we are going to formulate below. For this reason we consider the sets V m (Q [t , +00)) 
and M m (Q [to, +00)) . The former set is the set of column vectors of the dimension m 
and the last one is the set of square matrices of the order m G N. The components 
of the vectors from V m (Q [t , +00)) and matrices from M m (Q [t , +00)) are real valued 
piecewise continuous functions defined on the interval [to, +00) . We will also refer to them 
as to the linear spaces V m (Q [t , +00)) and M m (Q [t , +00)) , keeping in mind the sets of 
vectors/matrices mentioned above together with the addition operation and the operation 
of multiplication of the vector/matrix by a real number (see, for example, [161 P- 13-15]). 
Let (it, it) denotes the scalar product of the vectors it, it from V m (R) . We also define 

|^(t)||= SUp y/(lT (t),^(t)), 

te[a,b) ^32) 
(t) G V m (Q [t , +00)), 



\lt 



0,[a,b) 



sup 

t£[a,b) 



\\A (t)||| 



0,[a,b) 



sup 

t£[a,b) 



sup 



\A(t)lt\ 



r^ii / (33) 

\/A (t) G M m (Q [t , +00)) , ^ [a, b) C [t , +00) . 
It is easy to verify that the function |||-||| r a b ^ , defined by formulas (32), (33) is the 
norm of the linear spaces V m (Q [a, b)) and M m (Q [a, b)) . For convenience we will use the 
abbreviation 

IHIo = f IIHIIo,M) » when h b ) = [*o, +00) . (34) 
Let us fix a grid ^ on the interval [to, +00) and consider a set (a linear space) 
V m (Qlj [t Q , +00)) , which consists of the column vectors of dimension m. The components 
of the vectors from V m (Q* [to, +00)) are continuous real valued functions defined on the 
interval [t , +00) and continuously differentiable on each subinterval (tj_i,tj) , Wi G N, (we 
admit that there exist finite left and right derivatives at the points t i} i G N which can be 
unequal). Then we can define 

d 



I ^ (*)lli 



dt 



It {t) 



V¥ (t) G V m (Ql [t , +00)) 



(35) 



jit (t)||| = max {|||^ (t)\l , jit , Vlt (t) G V m (Qi [t , +00)) . (36) 

It is easy to verify that the functional ||j - [|| x (35) satisfies the conditions of seminorm in the 
space V m (Q* [t , +00)) , thus, the functional |||-||| (36) is the norm in V m (Q* [t , +00)) . 
We define the map P^: V m (Q^ [to, +00)) —> V m (Q [t , +00)) in the following way 

P h (It (t)) = it {U-x) , Vt G [*<_!, ti) , Vt G N, Vlt (t) G (Qi [t , +00)) . (37) 



Now we need to prove some auxiliary statements concerned with the operator P^,. 
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Lemma 2. Let for some &i > 0, u, v G Ei we have 

(u + 9iv, u + 6-lv) > (u, u) 
then for all 9', 6" such that 6' > 6" > ^ £/je inequality 



(u + 0'v, u + 0'v) > (u + 0"v, u + 0"v) . 



ZioZds trite. 



Proof. From (38) we would get 

(u, u) + 20i 3? ((u, v)) + 9\ (v, v) > (u, u) . 
From the last inequality, taking into account that Q\ > 0, we can obtain 

K«u,v»>-|(v,v>. 
Let us consider the first derivative of the function 

F(9) = (u + 0v,u + 0v) 

with respect to : 



(3? 



(39) 



(40) 



F' (0) = (u + 0v, u + 0v) e = 29£ ((v, u + 0v)) = 23? ((v, u)) + 20 (v, v) (41) 



We are going to prove that V > ^0i the expression (41) is nonnegative. Indeed, using 



the inequality (40) from (41) we would get 



F' (0) > -0i (v, v) + 20 (v, v) = (20 - 0i) (v, v) > 0. 

The last inequality means that on the interval [|0i,+oo) the function F (0) is nonde- 
creasing and V 0', 0": 0' > 0" > \6 l the inequality F (0') > F (0") holds true. The lemma 
is proved. ■ 

Lemma 3. For all u(t) G V m (Q^ [to, +°°)) , f or a ^ #2 such that Q\ > 02 and for the 



operator P h defined in (37) the following inequality 



|||P, (it (t)) - 0i (P h (it (t)) - It (t))\\ > |||P, (It (t)) - 2 (P, (It (t)) - It (t))||| (42) 
holds true. 

Proof. Let us fix an arbitrary n G N. As a first step we are going to prove the validity 
of the following inequality 

|||P, (It (t)) - 0i (P, (it (t)) - It (t))\\ 0>lt0;tn) > 

(43) 

> |||P, (It (t)) - 2 (P, (It (t)) - It (0)|«o l[tOjt „) • 
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Let t* G [to,t n ] be a point such that 

A = |||P, (it (t)) - 6 2 (P h (it (t)) - It W ;ino,[ to ,M 

= ||p, (it (t*)) - e 2 (p ft (r)) - ^ (t*)) || . 

If t* = t n then 

A=\\it (t n )\\ < |p, (it (t)) - e 1 (P h (it (*)) - (t))||| 0)Mn) v ^ g 



and inequality (43) is proved. 

Now we consider the case when t* ^t n . Without loss of generality we can assume that 
t* G Mi)- 

It is easy to see that A > \\lt (to)\\ , thus 



a = hp, (it (t*)) - e 2 (p h (it (**)) - it (t*))\\ = 

= \\!t(t )-e 2 (!t(t )-lt(t*))\\>\\!t(t )\\. 



(44) 



Inequality (44) means that the assumptions of lemma |2| would be satisfied if we put 
Ei = V m (Qlj [to, +oo)) , u = u (to) an d v = it (to) — u (t*) . Thus lemma [2] yields 

A=\\!t(t )-e 2 (lt(t )-lt(t*))\\ < \\lt(t )-9 1 (lt(t )-lt(t*))\\ < 

< \\P h (lt(t))-e 1 (P h (lt(t))-lt 

< \\P h (lt(t))-e 1 (P h (lt(t)) 



It 



0,[t ,*i) 
0,[to,t„) 



< 



So, inequality (43) has been proved. 

Due to the freedom of choice of the constant n G N inequality (43) yields (42). It 
concludes the proof. ■ 

Lemma 4. Let N (t, it) G M m ([to, +oo) x IR m )j and there exists a scalar function 

N (u) G C n (R) such that V (t, it) G [t , +oo) x IR m the following inequalities hold true 



k\ 



hi + . . .-\-krn — 

fci£NU{0} 



ki\ . . . k r , 



d k N (t,lt) 



du kl . . . du 1 ^ 



< 



d k 
du k 



N(u) 



Wk G 0,n, 



(45) 



then 



}i=0 



<A fc (iV(.);[|||^(t)|||o] 



i=0 / ' 



A k (N (t,-);[lt t (t)] 
Vlti (t) G V m (Q [t , +oo)) , zGOTfc, VA; G 0, n - 1 



(46) 
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A k (N (t, •) ; [Ji (t)]t ) - A k (N (t, ■) ; \P h (it, (t))) 
<hA k fiVW(-) x (•); fll^o (t) 1,..., Ift («)|| 



i=0 



< 



(47) 



V^i (t) G V m (Qi [to, +00)) , i G 0, k, VA; G 0, n - 2. 



Proof. Let us show that the assertion of lemma H] follows from lemma Q] To achieve 
that we need to put E x = (V m (Q [t , +00)) ; |||-||| ) , E 2 = (M m (Q [t Q , +00)) ; ||| 

■|lo) ' ^ n ~ 



Vm (Q [to, +°°)) ? £1 — (Vm (Qi [to, +00)) ; I -I) . It is easy to see that the operator P, (37) 
satisfies the following inequality 

(48) 



1^ (t) - P, (it (t))\\ < h \\lt C*) Hi , (t) G £ 1> 
moreover the inequality 

|p fe (it (t)) - e x (p h (^ (t)) - it (t))\\ > |||p, (it (*)) - # 2 (p h (it (t)) - it 



• 



v^i, e 2 eR:e 1 >e 2 >o, (t) g e x 



(49) 



is justified by lemma |3j 



In the considered partial case inequalities (48), (49) are equivalent to inequalities (23) 



and (25) respectively, and additionally the condition similar to condition (24) is triv- 



ially satisfied. Let us show that inequalities (45) implies (22). Indeed, for arbitrary 



it (t) , iti (t) G Ei, % G 1, k, Vt G [t , +00) we will have (see. PH p. 264, p. 279]) 

N§ } M(t))(^i (t),lt 2 (t),...,lt k (t)) 



mm m d k N (t, it) 

E U l,h (t) E U 2,i2 (*)••■ E U Ki k (t) 

i k =l W'i 1 ^"'i2 



! , _. I .J J ■, . i k dlli-, <)ll i , . . . fV// 



< E \ui, h (t)\... E \u k ,i h (t)\ 



ii=l 



a fe N (t, ^) 



m m 



<||^i(t)||...||^*(t)|| E ••• E 



du h du i2 . . . du ik 
d k N (t, it) 



t=it(t) 
< 



< 



(50) 



ii=l «i.=l 



mi^ 

i=i 



E 



du h du i2 . . . du ik 
d k N (t, it) 



k \ -\- . . . -\- k m — k 
fci6N(J{0} 



"a • • • • "771 ■ 



<9uf . . . du 



it=it{t) 



Tt=tt(t) 



< 
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< ( n ii^ 



k 

n 

i=l 



< n i^wiii 



iV(w) 



< 



«=IITt(t)l 



VA; G 0, n, 



ii=|||lt(t)| 



where iti (t) = [u^i (t) , . . . , u.j jTO (t)] T , Vi G 1, fc. Using definition]^ of the norm of fc-linear 



operator, from (50) it is easy to obtain inequalities (22). 



The presented above means that each condition of lemma [T] is satisfied, consequently 



its assertion, taking into account (35) and (36), implies inequalities (46) and (47). The 
proof is complete. ■ 

Corollary 1. Let N (t, it) G M m [c^'^ {[to, +oo) x IR m )j and there exists a scalar func- 
tion N (u) G C°° (R) such that V(t, u) G [to? +oo) x IR" 1 the following inequalities hold 
true 

d k N (t, it) 



E 



k\-\- . . .-\-forn — K 



k\ 

k\\ . . . k n 



du kl . . . du^ 



Then 



A k (n (t,-);[^ (*)] 



k 

i=0 



d k ~ 



< A k N 



, VA;GN|J{0}. (51) 



II It- 



I fc , 
U=o / ' 



(52) 



Vt G [t , +oo) , Vvti (t) G V m (Q [t , +oo)) , % G 0, fc, VA; G N (J {0} 



A, (N (t, •) ; [tf, (VOlo) " 4b (N (t, •) ; [It, (t)]t ) ||| o < 

< hA n (•) x (•) ; [ll^llio) , Vlt (t) G Kn (Qi [to, +oo)) 

Vt G fe-i.tj) , j G N, i G OJc, VA; G N|J{0} . 
The following lemma clarifies and generalizes lemma 2.2 from |12j . 



(53) 



Lemma 5. For an arbitrary scalar function N (u) G C°° (IR) , Vuj e 1, i 6 0,n-l the 
following equalities hold true 



1 d n |~/oo 



A n (JV(-) ;uo,ui,...,itn_i,0 
diV (w 



n! c?r r 



U=UQ 



i=0 



(54) 



Vn G N. 



T=0 



The convergence result for the FD-method applied to the Cauchy problem 
on the infinite interval. We will use the notation J (t, u) to describe the Jacobian 
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matrix of a vector- valued function N (t, it) it, N (t, it) G M m ([to, +00) 

dependent upon it (here we consider t as a parameter), i.e. 



x 



J(t,^) = N(t,lt) + 



(t,H 



It 



(55) 



[mi,m 2 , • • • ,u m } T . 

The following theorem gives the sufficient conditions for the convergence of FD-method 
applied to the Cauchy problem ([!]) on the interval [t , +00) in the sense of definition [TJ 
In particular, this means that if we are going to apply the FD-method to the Cauchy 
problem Q, which satisfies the conditions of the theorem stated below, then for an 
arbitrary interval [t , t + H) , H > we can use the uniform grid cu (JiJ) with the step h 
independent on H. It is easy to see that in such case the CPU time will be dependent 
on H only linearly (we assume that the computational complexity of the expressions for 
N (t, it) and 4> (t) has the same order on [to, +00)). The last fact is especially important 
when we are looking for the solution of the Cauchy problem ([!]) on a rather large interval. 

Theorem 2. Let the Cauchy problem Q satisfies the following conditions 



l)N(t,lt) 



E E 

p=0 il+...+in 



u\\..u^N h ... im (t), V{t,lt) E [t ,+oo)xR m , where 



N 



l\...l r , 



(t) G M m (C [t , +00)) , ik G N|J{0}, Vfc G l,m., and in addition 

to that there exists a sequence of nonnegative real numbers {-Bj}°^ such that 

00 

E UN 



"ll.-.lr, 



l! + ...+l m =p 

all mGI; 



(^)llo — Bp, G N|J{0}, and the series ^2u p B p converges for 

p=0 



2) ~${t) G V m (C%,+oo)), ~}{t) 



< K < +OO: 



3) there exists a constant a: < a G 1R such that \flt , it G V m (R) , Vt G [to, +00) the 
inequality 

(it, 3 (t, it) it) < -a (it, it) 

is valid. 

Then for an arbitrary initial condition lt G V m (M.) the solution of the Cauchy problem Q 
exists and is unique on [to, +00) . The FD-method for the Cauchy problem converges 
to the exact solution of the problem with the following error estimates 

C (h/R) p+1 



lt(t)-lt (t) 



lt(t)-lt (t) 



< 



< 



(p+l) 1+£ 1-h/R 
1 



h<R, 



R. 



(56) 



(57) 



i^i , + 1 ) 1+£ ' 

where the positive real constants C, R, e depend on the input data of the problem only 
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To prove theorem [2] we need several auxiliary statements, which will be stated below. 

of the- 



2) 



and 



3) 



Lemma 6. Let N (t, it) G M m {c t '^ ([t , +oo) x IR m )j and conditions 
orem [2] hold true, then for an arbitrary constant E\ > there exists a positive constant 
h = h (ei) G R such that for every grid u (jl]) wift h < h the solution (t) of the base 
problem ^ satisfies the inequality 

lll^ (0) (*)||lo<A*, 

(58) 

fi = max{||^o|| , ^} +ei. 

Proof. Let the assumptions of the lemma are satisfied. Let us set % = 1, fix an 
arbitrary t G [t , t\\ and multiply both sides of the system ^ on lt^ (t) . This results in 
the following (see pH p. 284 - 288]) 



= (it® (t) , (N (t, it® (t)) ^(°) (t) - N (t, ~$)~$)) + 
+ <^(°) (t) , (N (t, ^(°) (t )) - N (t, (t))) ^ (t)> + 



(59) 



+ ( t*(°> (f) ,"£(*)) = <^(°) (t) , J (t, tfy (t)> + 



+ (t - to) /^ (0) (t) , T (t, lt 2 , it® (t)) Qtf <°) (t) 



where lt x = erf® (t) ,lt 2 = it® (t ) + 6 2 (lt^ (t) - (*„)) ,6 k e[0,l], k = 1, 2, 
t G [to,t], 



T(t,^,^) 



dN(*,lt) , dN(t,lt) dN(t,lt), 
-v, — ~ v — v 



dui 



du r 



Taking the norms in the last equality of (59) we get the following 
1 d 



2dt 



|^(°) (tW < -a II ^(°> (t)|| 2 + k (t)|| + 



+h 1 \\it<n (t)\\\\T (t,it 2 ,^ (0)11 x 

x {pM (t) || ||N (t, it (t )) || + k} , h = h - 1 . 



(60) 



Let us fix an arbitrary E\ > and, taking into account (58), put 



/ii = max 



{11*11, 



+ 



(61) 
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It is easy to verify that in this case a\i\ — k > 0. We are going to prove that if 

ajii — k 



0<ht<h<h 



72 (ji) (7i 0) V + «) 



(62) 



where 



71 (fi) = sup ||N(t, < +00, 

iteBft,, te[t ,+oo) 

12(1^) = __ sup ||T (*, "tf, < +00, 
itjlfeBp, te[t ,+oo) 



(63) 



= {it e V, 



then 



sup 

te[t ,ti 



"II — A*} 



|^ (0) (t)|| <m<ti. 



(64) 



Indeed, let us assume that under condition (|62|) inequality (|64|) is not valid. So, there 

(65) 



exists the real numbers £*, t* : t < t* < t* < t\ such that 

m < \\it (0) (t)\\ < /i, vt g (t*,t*] , p (0) (t* 



Then from (60) and (62) Vt G [t*,t*] we obtain 



1 d 
2dl 



|^ (0) (t)\\ 2 < ||^ (0) (t)||(-a//i + K+(a/ii-K)) = 0. 



(66) 



Inequality (66) shows that the function [["u^ ) (t) || is nonincreasing on [t*,t*], hence 
|| u ( -°- ) (t) || < fix, Vt G [t*,t*]. The last inequality contradicts with assumption (65). 
Thereby, inequality (64) is valid. 

Similarly, if we consider the base problem ^ with i — 2, it is easy to prove, using 



(64), that under condition (62) the inequality 

|^ (0) (t)\\ <fii<H 



sup 

te[h,t 2 ] 



holds true. In general, we can prove that condition (62) together with the assumption 
\\lt^ (ti-i)|| < /ii implies the inequality 



sup 

*e[ti-i,t 4 ] 



|^ (0) (t)\\ <fii<fM. 



Thus, using the principle of mathematical induction, we obtain 



sup 

te[to,+oo] 



P {0) (t)\\ <fii<fi. 



(67) 



From (67) we obtain inequality (58) for the h defined in (62). The lemma is proved. 
Lemma 7. Let J (t) be the matrix from M m (C [t±, tz]) , t\ < tz and 



(it 



J (t) it) < -a {it, it) , < a G E, \/lt G V m (R) , Vt G [ti, t 2 ] , 



(68) 
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then the matrix Q (t) G M m (C 1 [ti,t 2 ]) which is the solution of the Cauchy problem 
d 



dt 



n(t)-j(t)n(t) = o, n(h) = E, te[t 1 ,t 2 \ 



(69) 



satisfies the following inequalities 

\\Q{t)\\ <e- a ^\ te[t u t 2 ], (70) 
HO" 1 (t)|| < exp ((* - ^) HI^C^IIIo,^^)) , te[t h t 2 }. (71) 
Proof. Let the assumptions of the lemma hold true and Q (t) denotes the solution of 



problem (69), then V¥ G V m (R) we have 

1 d 



1 d 



, )<h ^(t)^t = ~^(t)^^(i) it) 
(J (t) Q (t) it, n (t) it) < -a (Q (t) it, tt (t) it) = -a \\Q (t) lt\ 



(72) 



From (72) we obtain 



dt 



\n(t)H 



— > II 2 



\n(t)H 



— > II 2 



< -2a. 



(73) 



Integrating the both sides of inequality (73) from ti to an arbitrary fixed point t G (t\,t 2 \ 
we get the inequality 



In (t) lt\f\ - In (\\lt\f\ < -2a (t - h) . 
So far as inequality (74) is valid for any arbitrary vector it, we have 



(74) 



sup In = r— 2 — J = In sup , 2 

ltey m (K) \ ||«|| / wev m (R) II u I 1 " 

|2\ 



(75) 



In (I|S^ (*)lr) <-2a(t-t!) 



Inequality (75) implies (70). 



Now we are going to prove inequality ( |7T| ). To do that we need to differentiate both 
parts of the identity Q (t) Q^ 1 (t) — E,\/t G [t 1 , t 2 ] with respect to t. Then we get 



d 
dt 



rr 1 (t) = -fi- 1 (t) j (t) , n- 1 (to = e. 



(76) 



Integrating the differential equation (76) from t\ to an arbitrary point t G (ii, £2] we obtain 



t 



(77) 



fl 
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Taking the norms in the last equality we come to the following inequality 



Z I 

|n~ 1 (*)||<PII + f ||n- 1 (e)||||J(OI|de = i + f p- 1 (0|| 11.7(011 (re 



ti 



Applying the Gronwall-Bellman inequality (see [TU p. 108]) to (78) we obtain (71). The 
lemma is proved. ■ 

Proof (of the theorem [2]). Let the conditions of the theorem hold true. The 
existence and uniqueness of the solution of the Cauchy problem ^ on [t , +oo) was 
proved in [TU p. 286]. 

Let us fix an infinite grid cu Q. Then using a scalar parameter r G [0, 1] we consider 
the following generalization of problem ([I]) 



d 



it (t, r) - N (t, it (ti_i, r)) it (t, r) - 



-r {N (t, it (t, r)) - N (t, it (t»_ ls r))} it (t, r) = ^ (t) , 



(79) 



rG [0,1], te [U-uU), V* GM, 
it (t , r) = lt , [it (t, r)] t=ti = 0, zGN, Vr G [0, 1] . (80) 

It is easy to see that if r = 1 then problem ( 79 ) , ( 80 ) transforms to (JTJ . 

We assume that problem (79), (80) has a unique solution u (r, t) , which can be ex- 
pressed in the series form 



it (r, t) = J2 ^ W (*) » Vr e [0, 1] , Vt G [t , +oo) 



(81) 



i=0 



where «W (t) G V™ (Q* [to, +oo)) , i G M|J {0} (see notations on p. 10), furthermore 



d 00 d 00 

It (r, t) = J2 t*-^® it) , Vr G [0, 1] , Vt G |J (t^, U) 



Of 



i=0 



<9 



lt(r,t) 



^2) 



^r^M-^fi), VrG [0,1] Vj'GN|J{0}. 



t=ti i=0 



Taking into account assumptions (81 ), (82) and putting r = in (79), (80), we obtain the 



base problem (|5| to define unknown term it'- -' (t) . Similarly, if we substitute representa- 



tion (81) into (79), differentiate the obtained equality with respect to r j times (j G N) 
then divide both sides of the equality by j\ and finally put r = 0, we obtain the recursive 
system of Cauchy problems (|6| to define unknown terms u^) (£) ) \fj g pj. 
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Let us prove that the FD-method for the Cauchy problem ([!]) converges in the sense 
of definition [T] To do this we have to rewrite ^ in the equivalent form, namely: 



d 
Jt 



^0) {t ) - (N (t, (^0) + T, (t, (««_!))) ^ (t) 



53) 



= T, (t, t^W (t)) If W (t^) - T, (t, ^(°) (t^)) ^°' } (t) + F&) (t) 
t G [ti-i,ti) , z, j G N, where, V^T" G (M) 



It 



34) 



Taking into account notations (55), we can rewrite equation (83) in the following form 

* d 



d 
Jt 



^ (t) - 3 t (t) {t) = T t (t, J -^(°) (0 d£) {t^x) - 
Ji (t) = J (t, lt<® (Vi)) , *e i.jeN. 



55) 



Then if we add to equation (85) the initial and matching conditions 

^')(t ) = ^o, P (i) (*)], =t . = o, yeN, 



56) 



we would get the recursive system of the Cauchy problems with respect to the unknown 
vector-functions it) , which is equivalent to system (|6]). 

Let us fix an arbitrary j G N. Then the solution u^) (£) Q f the Cauchy problem (85), 
(86) can be expressed in the form 

it® (t) : 



ti-1 



e d 



e d 



k-i U-i ar l 



(87) 



where the matrices Qi (t) G M m (C 1 [ii_i,£t]) are the solutions of the Cauchy problems 
d 



dt 



■Qi (t) - Ji (t) a (t) = 0, Qi (t<_i) = £, iGN, 



^i(t,0 = fii (*) n; 1 (o. 
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Lets use lemma [6] and assume that for some fixed constant E\ > 0, a maximum step h 
of the grid u) Q satisfies the inequality h < h, where h = h(e\) is a constant mentioned 
in lemma [6j Since that, accordingly to the assertion of lemma [6j we can assume that for 
the vector-function lt^ (t) that is the solution of the base problem ^ inequality (58) 
is valid. Whereas the vector function (t) depends on the grid Q (0) it is worth to 



emphasize that the right side of estimate (58) does not depend on u). After that using the 
following notations 

dN (t, it) 



B = max < ^2 

\\~u <Ai I i=l 



dui 



C = /i max < ^2 



N (t, it) it + (f> it) 
dN (t, it 



< +oo, 



dm 



< +oo 



and taking into account the evident inequality 
\\T t (t,lt(t))^ 



™ <9N (t, it 

?>w-sr- 



— > 



i=l 



dN (t, it) 




dui 


•tf=-*e>)(ti_i) 



t?=tft°)(ti_l). 

^(*)loll^(*)llo. V?eN - 



< 



(t) , (t) e V m (C [t , +oo)) from (|87| we would obtain 



< 



1 + 



i=i 



dN (t, it) 



hi m 
^ 8=1 



dui 
dN (t, it) 



dt 



!tW (t) 



l^ (0) W|lo 

lZ > =^(0)(t l _ 1 ) 

Will < 



lt^ (t) 



+ 



0,[*<_x,ti) 



< 



tiB 



\lt^ (ti-t) 



dt 



lt^ (t) 



+ h t \\\F^ (t)||| o 



0,[t 4 -iA) 



Similarly, from (87) we obtain 



\!tW (ti)\\ < 



e —hiCe _|_ 



ti?B 



\lt {3) (U-JW+hi \\\F^ (t)\\\ Q + 



dt 



it® (t) 



i e N. 



59) 



(90) 
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On the other hand, from ^ we find 

tf <fl (t) = U t (t) It® (^0 + 17, (t) / f/r 1 (0 F U) (0 + 

+^ (*) / f/r 1 (0 t, (e, ^ (0) (0) fe-i) . 

where the matrix C/j (t) G M m (C 1 [t^_i, t,;]) is the solution of the Cauchy problem 



-Ut (t) - N (t, ^ (0) (*<_!)) (t) C/, (t) = 0, C/i = £7. 
If we differentiate with respect to t both sides of equality (91) we would obtain 
d 



dt 



^ (t) = N t (t) Ui (t) + T t (t, ^(°) (t)) + 



+ iv, ( t ) ^ (t) j u,- 1 (o t, (e, (o) # 

ti—i 

+N t (t) u (t) j ur 1 (0 (0 d£ + F«) (t) , t e [U-t, U],ie N, 

ti— 1 

where iV, (t) = N (t, "tf (0) . From Q we get 

ffliV,(t)|| exp{/i|^(t)| } + 



< 



0,[ti_i,t<) 



m 

+ E 

i=i 







dui 





ll ^(t)||| + / i |||iV,(t)||| exp{2/ i |||iV i (t)||| }x 



8=1 



<9N (t, if 



ll^Willo 



tf® (ti-l)\\ + 



where 



+ (h IN, (t)|| exp {2/. |JV, (t)| } + 1) |||FW (t)||| < 

<p||^) (w)||+q|||fw (t)||| 



P = JVe^ + Cg, Q = hNe 2hN + 1, iV = max |||N (t, ^) ||| . 

*M — ^ 



From inequalities ( |89| ), (90) by virtue of (93) we obtain 

\^ {3) (t)\\\ , [ti _ uU] < [1 + (ti-i)H + /il> (t)| 



( fi )|| < [ e -«fc + C**— x)|| + KD \\\F^ (t)\\\ Q 
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where 



B + CP CQh 
E = , D = l + ——. 



Let us make in (96) the following substitution 



\ fU) wlllo> ' 



(97) 



then the system of inequalities (96) turns to the following one 



y . < [ e -«* + Etil] Vi _ x + htD, y = 0, i G N. (98) 

Now we are going to prove that for every sequence of sufficiently small constants hi > 
0, ieNa sequence {yi}^ which satisfies inequalities (98) is bounded. So far as for all 
hi > the following inequality holds true 



-ahi 



+ Ehf < 1 - ah, t + hn E + 



a 



1 - hi[ a - hA E + 



a' 



the assumption 



yields the inequality 



hi < 



a 



2E + a 2 ' 



i e N, 



(99) 



-ahi 



+ Ehi<l-ahi + ht[E+ — \ < 1 



hid 



Taking into account the last estimate it is easy to verify that under conditions (99) the 
recursive sequence 



Y; 



a 



l-h i -)Y i _ 1 + h i D, J6N,y„ = 0, 



(100) 



is dominant for every sequence {yi}^L which satisfies inequalities (98). If we make the 
substitution 

2 2D 



t)i~, Yi 
a 



a 



[1-Zi), iE N, 



in (100) we would obtain the following recursive sequence 

Z i = (l-f) i )Z^ l , ieN, 2-0 = 1. 



Thus 



and the assumptions 



Zi = Y[(i-f) p ), ieN, z = i, 
P =i 



V/i,- : < hi < h = min 



4 a _ 



a 2E + a 2 



(101) 
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implies the estimate 



< yi < Yi < -D. 

a 



(102) 



Using estimate (|102|) and taking into account (|97|) from inequalities (|93|), (|95|) we can get 

d w ., , i /7' 



eft 



^ (t) 



< 



a 



+ Q)\\\F^ (t 



\^ U) (t)\\\ < ( (l + Eh 2 ^) ^ + hDj \\\F^ (t)||| 



Let us put 



r4PL> 

(7 = max < h <5, ( 1 + Eh 



AD 



a 



hD 



then from (103) and (104) we obtain 

W||| <a\\\F ij) (t)||| . 



(103) 



(104) 



(105) 



In the following inequalities we use the assertions of corollary [T] and lemma [5] because 
it is easy to see that their assumptions are satisfied if we put 



N{u) = J2 B i u ^ Vn G 

j=0 



(106) 



Indeed, using the Leibniz formula for the n-th derivative of the product (see, for example 
[TTt p. 268-272]) we can check that the function N (t) (106) satisfies conditions (51): 



d k B p u p 



[B v vP 



,(*) 



u= \\~v 



u= 1 1 it || 

(ft) 



> E ll|N,,., m (t)||| ^ + -^ 

\ i 1 +...+i m =p 

E , , k \ , E |N, 1 .., m (t)||| (« 



> 



u=|M 



E 



jfc! 



fci+..!+fc m =fc hi... k m \ i 1+ .^i m=p ^ n '" %m ^ ^° du^du^ 2 ■ ■ ■ du^ 



d k (u^ul 2 ...u%) 



> 



k 1 + ...+k m =k r»l- ■ ■ ■ ""m- 



E 

ii+...+i m =p 



d k N n ... lm (t)u^...ui 



du kl du k2 ...du^ 
V If = [vi, v 2 , . . ■ , v m ] T e V m (R) . 
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From (105) we have 



f j ~ 1 



(t)||| < at e Vp (•) ; WlllE) IIP (P) Will + 



E (•) >< (•) ; [|||^ (i) (*)IK: 1 ~ P ) IP 00 Will + 

p=0 v 7 

III ^ (0) Will dP 



N ( E ri W| 

,i=0 



(107) 



(t)\ 



r=0 



Let us put 



(t)||| = Vj , jeN. 



(10* 



in (107). Then, changing Vj by Vj and the inequality sign to the equality one we obtain 
the following system of equations with respect to Vj 



V S = <t\ E Aj-p (N (•) ; [VrJ) V p + E Aj-i-p (n> (■) x (■) ; [V] 



P =i 



+ 



V d j ( ~ / 



'ins 



p=0 



T=0 



i-i-p 

8=0 



where Vq = /x, or 



(7 



1 + O-VniV' 



(Vo)1p 



3-1 



(109) 



+ E Aj+ P (N> (■) x (•) ; m^ P ) V p },jeK 

p=0 



The sequence , which satisfies the system of recursive equations (109), is dominant 

for the sequence . Using the generating functions method, we are going to prove 

that for the h sufficiently small the assumptions (81), (82) hold and the power series 

oo 

z % Vi possess a nonzero convergent radius. 
i=0 



From (109) we have 

g(z)-v 



1 + aV N' (Vq) 



g(z)(N(g(z))-N(V )) + 



(110) 



)N'(g(z))} 



where 



+zg [z 



3=0 



(111) 
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To prove that the function g (z) , which satisfies equality (110), exists or, in other 
words, has nonempty domain which includes some neighbourhood of the point z = 0, we 



have to express z from (110) by means of g: 

1 



z{g) 



9 2 N' (g) 

V <g, S 



E(g-V )+[N(V Q )-N(g))g 



(112) 



1 



a 



+ V N' (V 



It is easy to see that the function z (g) ( |112 ) is defined and continuously differentiable in 
some neighbourhood of the point g — Vq. To prove that there exists an inverse function 
g = g (z) defined in some neighbourhood of the point z = 0, it is sufficient to check that 



z' (Vo) > 0. The last inequality follows from the expression for z{g) (112) 

z{g)-z(V ) 



z'(V ) 



lim — = 

g^V g 2 N t ( g ) 



lim 

g-*Vo g — V 

N(g)-N(V ) 
9-V 



(113) 



crV 2 N> (V 



> 0. 



It follows from inequality (113) (see [HI p. 87]) that there exists an inverse function 
/ (z) which is holomorphic in some open neighbourhood of the point z = : C C C. 
Assume that -Ri > is the largest constant for which the open ball = {zeC: \z\ < 



Ri} satisfies inclusion C C. Now we are going to prove that series (111) converges 
at the point z — R\. Suppose that this is not true and series (111) diverges at the point 
z = R\. Then we have 

lim f(z) = +oo. 



But, on the other hand, for all z from (— Ri) equality ( 110 ) holds. And we immediately 
get the following contradiction 



l = f(z)-v + 
(l + v ) z [\\ g '\\ 0il f( z ) + N(f(z)) + N'(v )v -N(v )] 

+ 7^ = ► +00, 

fyz) -v 

as z — > R\ — 0. This contradiction implies the inequality f(Ri) < +oo and 

C 

3 ~ U + 1) + 

for some positive constants s and C, which do not depend on h. 



(114) 
(115) 



If we return to notations (108) we would obtain an estimate 

c ( hy . . 



will < 



(j + iy 



(116) 



2G 



and the following sufficient condition for series (fff ) be convergent 

h 

ft* 1 ' 

So, we have proved that for the sufficiently small maximum step h > of the grid cD 
Q, namely, 

h < R = min {X-Ri} (117) 

(see (101)) assumptions (81), ( [82] ) hold true. Hence, in accordance with definition [IJ the 
FD-method for the Cauchy problem Q converges. Now we have to prove that it converges 
to the exact solution of problem Q by the same definition. 

Let us assume that the grig u Q satisfies condition (117). We define (see notations 
on page 



10) 



It easy to see that 



-&(t) = Y, ~J U) (*) e Vm (Qi K +00)) 
3=0 



d 00 d 



(118) 



dt 



j=0 



Let us prove that vector-function (118) is the solution of the Cauchy problem ([I]). For this 
reason we need to sum the equation of the base problem ^ with equations ^ Vj G N. 
It results in the following equalities 



d ^(t) -EE (N (t, •) ; p« (f)]£) « (t) = ^ (t) 



^(HO)-EE (N ft, •) ; [tf M ft)] j^j) « ft) = (*0 
j=0 p=0 v 7 



;ii9) 



VzgnuW, vtG u (t,-!,^)- 

i=l 

Using the theorem about the substitution of the power series into the power series (see, 
for example, [201 P- 485]), it is easy to obtain the following equality 



(OO \ oo 

*,5V^ W (*))£ T<1 * W (t) 
i=0 ' i=0 



i=0 p=0 
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Wt G [0, 1] , Vt G [t , +oo) . Taking it into consideration, we can rewrite equalities (119) in 
the following form 



-IT (t) — N(t,v (t)) v (t) — <p (t) ,t E u fe, t f _o , 
at i=1 

^ (*i + 0) - N 17 (*,)) 17 ft) = ~$ (U) , i G N U {0} • 



(120) 



The uniqueness of the solution u (t) of the Cauchy problem Q together with equalities 
(120) imply that as far as 17 = if the identity lt(t) = lt (t) , Vt G [£;_!, , 
i G N holds true. Thus, the evident equality v (t ) = it (t ) = u implies the identity 
17 (t) = if (t) , Vt G [t , +oo) . 

Combining ( 116[ ) with (117) it is easy to obtain estimates (56), (57). The theorem is 
proved. ■ 

Numerical example. Let us consider the following Cauchy problem 



d 
It 



u x (t) + (u\ (t) + u\ (t) + 1) u x (t) - Ul (t) u 2 (t) 



2 sin (t) - -sin (2t) + cos (t) , 



d 



u 2 (t) + (uj (t) + u\ (t) + 1) u 2 (t) - Ui (t) u 2 (t) 



2 cos (t) - -sin(2t) - sin (t) . 



w 1 (0) = 0, m 2 (0) = 1, tG[0,+oo). 
Using the notations of system (JTJ) we would have (it = [ui,u 2 ] T ) 



N(lt) 



U 2 — u\ — «2 — 1 

u\ — u\ — u\ — \ 



1 















1 







1 







1 









+ Ml 






+ u 2 






+ u\ 






+ u\ 









1 





1 











1 







1 



(*) 



N (u) = l + 2u + 2u 2 



2 sin (t) - -sin (2b) + cos (t) 
1 

2 cos (t) - -sin(2i) - sin (t) 



It is easy to verify that the Jacobian matrix of the vector-function 

VI (it) = N (it) it = [ui (u 2 - u\ - u\ - l) ,«2 (ui -u\-u\- 1)] : 



;i2i) 



(122) 



(123) 



(124) 



(125) 
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has the form 



Jilt) 



u 2 — u 2 — 3u\ — 1 
u 2 - 2u 1 u 2 



U\ — 1U\U 2 



U\ — 11 



3u 2 2 



1 



Applying the Sylvester's criterion to the symmetric matrix J s (it) = - ( J (it) + J T (it)) , 

we obtain that the matrix J s (it) is negative defined \/lt G V 2 (R). If we find the extremum 

Tr (J s (it)) 

values of the function - — — r— , as the function of two variables u±,u 2 , we would get 

det (J s (H)) 

the inequality 



Tr(J s (lt)) Ai (H) + A 2 (IT) , 

> ; — = — ; > -2.15, VH E V 2 



det (J s (It)) ~ A x (lt)\ 2 (lt) 



(126) 



where Ai ( u ) , A 2 (¥) < are the eigenvalues of the matrix J s (u) . From inequality ( |126[ ) 
it follows that 

min{|Ai(^)|,|A 2 (^)|} >a = (2.15)" 1 « 0.47. (127) 



Similarly, if we find the extremum values of the function 
the estimate 

||0(f)ffl o <« = 2.9. 



we would come to 



:i28i 



So, as it was showed above, the Cauchy problem (121), (122) satisfies the conditions 



of theorem |2| with the constants a (127) and k (128). It is easy to see that the exact 
solution it* (t) of the problem is 

it* (t) = [u\ (t) , u* 2 (t)] , u\ (t) = sin (t) , u* 2 (t) = cos (t) . (129) 

As a first step we try to use the ADM to find the approximate solution of problem 



( pi) , ( p2| ) (see, for example, [I]). 

For this reason we need to express the vector-function 9T (it) ( 125 ) as a sum of a linear 
component and an essentially nonlinear one: rfl(lt) = m 1 (it) - E^lt. We looking for 
the p-th approximation of the exact solution of the Cauchy problem (121), (122) in the 

form of the partial sum u a if) = [uai (t) , UA2 (t)] = ^2 It a f) where unknown vector- 

i=0 

functions lt^ (t) are the solutions of the following Cauchy problems 

d If (t) = -E^ltf (t) + -$ (0 , Itf (0) = [0, If , 

(130) 



dt 



d 
dt 



Itf (t) = -E^ltf (t) + ~ (rVli (fW? (0)) 



r=0 



-E^itf (t)+A-x (W),^? (0 



i-1 
k=0 



■^w (o) = 0, i e N. 
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b) 



Figure 1: The Cauchy problem dm}, ((ml), ADM. □: (t) , <>: u Al (t) , +: u Al (t) , o: Uj4l (t) , the 



solid line represents the exact solution (129), a) i = 1, b) i = 2 



The results shown on figure [I] tell us that the ADM ( 130 ) is divergent on the segment 
[0, 2] . 



Applying the FD-method to the Cauchy problem ( 121 ), ( 122 ) we used the uniform grid 

with the step h = 1/5. We have found five approximations u (t) , p G 0, 4 of the exact 
solution on the segment [0, 6] . Also we have calculated the absolute error values 



5f (t) = u t (t) - Ul (t) 



6 {p) (t) 



lt(t)-^ (t) 



1,2, pe0,4, 



p e 0,4. 



The numerical results showed on figure [2] confirm that the FD-method for the Cauchy 



problem (121), (122) converges with the exponential rate to the exact solution of the 
problem. 
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